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ABSTRACT 

We present the first results from the Distant Radio Galaxies Optically Non-detected in the SDSS (DRaGONS) 
Survey. Using a novel selection technique for identifying high redshift radio galaxy (HzRG) candidates, a large 
sample is compiled using bright (5i.4ghz > lOOmJy) radio sources from the Faint Images of the Radio 
Sky at Twenty centimeters (FIRST) survey having no optical counterpart in the Sloan Digital Sky Survey 
(SDSS). Near-infrared (NIR) JC-band imaging with the FLAMINGOS instrument on the 4-meter telescope at 
Kitt Peak for 96 such candidates allows preliminary identification of HzRG candidates through the well-known 
K — z relation, and these objects will subsequently be observed spectroscopically to confirm their redshifts. 
Of the initial candidates, we identify 70 with magnitudes brighter than K « 19.5, and compute limiting 
magnitudes for the remainder Assigning redshifts based on a linear fit to the K — z Hubble diagram gives a 
mean redshift for our sample of z = 2.5 and a median redshift of z = 2.0, showing that this method should be 
very efficient at identifying a large number of HzRGs. This selection is also sensitive to a previously unseen 
population of anomalously red radio galaxies (r ~ K > 6.5 — 7), which may indicate significant obscuration 
at moderate redshifts. These obscured objects can be used to test the completeness of QSO surveys to the 
effects of reddenning. More than ten percent of our sample falls into this category, which may represent a 
sizable radio loud population missing from current optically selected AGN samples. We additionally identify 
479 bright Extremely Red Objects (EROs) in the fields surrounding our HzRG candidates, to a magnitude 
of K = 17.5, within a non-contiguous area of 2.38 square degrees. This constitutes a small overdensity of 
EROs surrounding the radio galaxy candidates over random fields, and we see possible evidence for a physical 
association of the EROs with the radio galaxies. However, examining the clustering of all K < 19.0 galaxies 
around the radio targets reveals no evidence of a global galaxy excess, strengthening our conclusion that the 
EROs trace an overdensity not evident in the overall galaxy population. 
Subject headings: galaxies: evolution — radio continuum: galaxies 



1. INTRODUCTION 

Radio galaxies ha ve long been used to probe the epoch 
of galaxy fo rmation (Lillv et al. 1984; Graham & Dev 199^ 
iBlundell et al. 1998; Stern et al. 1999; Jai'vis et al. 2001). Lu- 
minous radio galaxies are known to be highly biased relative 
to the underlying dark matter, residing in the most overdense 
regions of the Universe. Under the standard Cold Dark Mat- 
ter hierarchy for galaxy formation these galaxies should be 
the first systems to collapse and, therefore, the site of some 
of the most evolved stellar populations. As the likely hosts 
of the earliest star formation, it may also be possible to probe 
the epoch of reionization (Barkana & Loeb 2005). Isolating 
a sample of radio galaxies at high redshift would allow us 
to probe the physical processes that drive the formation and 
evolution of structure and the timescales that govern star for- 
mation in the early universe, as well as to distinguish between 
hierarchical and "down-sizing" iCo wie et al. 19961 ) formation 
scenarios (e.g., Rocca-Volmerange et al. 2004; de Lucia et al.' 
|2p05), and also the role of feedback on the local radio galaxy 
environment ( Vardoulaki et al. 2005; de Lucia et al. 2005). 

While radio surveys can address many fundamental ques- 

' Based on observations collected at Kitt Peak National Observatory, Na- 
tional Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc. (AURA) under cooperative 
agreement with the National Science Foundation. 



tions in cosmology and galaxy formation the numbers of radio 
galaxies identified at high redshift remain small in comparison 
to low redshift surveys ( B rand et a ni2005tlMagiiochetti et alJ 
l2004h . Without large, statistically complete and homo- 
geneously selected samples we cannot hope to constrain 
hierarchical galaxy formation models without the concern 
that sample variance might bias our analyses. For ex- 
ample, at redshifts z > 3 extensive optical and near- 
infrared (NIR ^ campaigns have yielded less than one hun- 
dred galaxies (Van Bruegel et al. 1998: de Br euck et alJ200l|; 
Vardoulaki et al. 2005). The reason for the paucity of these 
samples comes from the necessity of surveying large volumes 
to identify the most massive systems. Given the broad redshift 
distribution of radio galaxies (e.g., Dunlop & Peacock 1990), 
large numbers of radio targets must be observed in order to 
extract the high redshift component. 

In this paper we describe a novel selection technique that 
utilizes existing optical and radio surveys to overcome these 
challenges in order to isolate high redshift candidates for fol- 
low up in the NIR. We show that this approach gives more 
than a factor of ten increase in efficiency over blind radio se- 
lected surveys. The resulting combination of NIR data with 
existing optical measurements from SDSS allows us to study 
the environment of the high redshift radio galaxies. Early in 
the development of infrared observations it was noted that 
EROs are often found in the vicinity of high redshift ob- 
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Fig. 1.— dN/dz based on the IDunlop & PeacocU IT99(1) model for 
•S'i.4GHz > 100 mJy radio sources. This shows the broad redshift distribu- 
tion expected for bright radio sources. The y-axis scale is in arbitrary units. 



jects ("McCarthy. Persson & West"1992>. rCiraham et aUflOOl 
iDev. Spinrad. & Dickinson 1995). More recently, targeted 
searches have b een undertaken to search for clustering of 
galaxies ('Hall et al. 2001) and extremely red objects (EROs) 
around high redshift quasars and radio galaxies (|Cimatti et al. 
l2000t IWoldetalJ l2003: Zheng etal. 2005, and references 
therein). The focus of these investigations has been whether 
the increased surface density of EROs is physically associ- 
ated with the radio galaxy or quasar, or whether it lies in the 
foreground, possibly as a cluster that may have gravitational 
lensing effects on the target. 

Also of interest in a near-infrared radio survey is the num- 
ber of obscured quasars. It is well know that optical selection 
is sensitive to a wide range of effects that bias samples against 
the detection of heavily obscured galaxy populations. Indeed, 
heavily obscured quasars may account for a significant frac- 
tion of the total popu lation (Webster et al. 1995; White et al. 
^003: Gli kman et al.l 120041) . Our selection criteria result in 
the detection of a population of bright {K < 17.5), red 
(r — K > 6.5) objects that may fall in to this category. 

The structure of this first in a series of papers is as follows: 
The target selection procedure is presented in §|2l A descrip- 
tion of the observations, astrometry and photometry is given 
in §|3] The results are presented and discussed in §0] and a 
summary and discussion of future work is presented in §|5] 

SDSS magnitudes are in the AB system, while we will 
use Vega magnitudes for all ii'-band measurements for easy 
comparison with the high redshift radio galaxy literature. 
A flat Lambda cosmology is assumed throughout, as in- 
dicated by numerous recent measurements, with Hn = 
70kms-iMpc-\ = 0.3, VIa = 0.7 (e.g.. lSpergel et al.1 
12003,) . 

2. TARGET SELECTION 



The main obstacle to identifying high redshift radio galax- 
ies is screening out the low redshift foreground. Figure ^ 
shows a dN/dz distribution for bright (5i.4ghz > 100 mJy) 
radio sources based on model radio luminosity functions and 
assum ing pure luminosity evolution ( Rowan-Robinson et al] 
ll993t [Dunlop & Peacock 1990). With a broad redshift peak 
al z ^ 2 extracting only the high-z galaxies through blind 
spectroscopic follow up of radio surveys is inefficient. Our 
goal is to eliminate the low redshift contamination through 
the inclusion of multi-wavelength information. We begin with 
the well known K — z Hubble relation for radio galaxie s 
dlJllv et alJl9 84 HVan Bniegel et alJl 99a lj arvis et alJ200lD . 
which shows the strong correlation of K-b&nd apparent mag- 
nitude and redshift. If we assume that this relation holds for 
all bright radio galaxies, then we can use model galaxy col- 
ors to predict the optical properties of these galaxies as they 
evolve. 

Figure|2]shows the r — K, g — K, and i — K color-redshift 
diagrams for two sets of galaxy models generated with the 
spectr al synthesis code PEGASE {Fioc & Rocca-Volmerange] 
The upper curves are color tracks for an instantaneous 
burst model at solar metallicity and assuming no extinction 
for formation redshifts of Zfoim = 10, 5, and 3, while the 
lower curves are for zero initial metallicity with an exponen- 
tially declining star-formation model with an e-folding time 
of r 0.5 Gyr and PEGASE "spheroid" type obscuration at 
these same formation redshifts. We choose solar metallicity 
for the instantaneous burst model to match the r — K color of 
an elliptical galaxy at low redshift, as the PEGASE models do 
not update the metallicity in an instantaneous burst. Metallic- 
ity does evolve in the models of ongoing star formation, hence 
the assumption of zero initial metallicity for the r = 0.5 Gyr 
e-folding star formation model. Although the spheroidal ob- 
scuration makes the ongoing star-formation model slightly 
redder than a typical elliptical galaxy at low redshift, we in- 
clude it to approximate the effects of dust during the starburst 
(i.e. PEGASE does not allow for the destruction of dust). We 
would expect this model to be accurate at high redshift, but 
too red at low redshift. Since the low redshift models are ex- 
cluded in both cases, this is not a concern. Shown for refer- 
ence on each of the color redshift plots is the line representing 
the color of a point source at the 5 a limiting magnitude (AB) 
of the SDSS filter (.gum = 23.3, = 23.1, im„ = 22.3) 
with the /f-band magnitude assumed from a linear fit to the 
K — z diagram, given by: K = 4.62 log(z) + 1 7 .2, de - 
rived from a fit to the galaxies in Van Bruegel et al. ( 1998"). 
Sources below the line on the diagram would be detected in 
each SDSS band, while sources above would not. Note that 
there will be some scatter in the cutoff due to the scatter in the 
K — z Hubble diagram. The exact nature of the objects pass- 
ing the magnitude cut will depend on the color dependence of 
the scatter in th K — z relation. Combining the information in 
the K — z diagram and Figure|2]we see that all of the optically 
bright galaxies are either at low redshift, or galaxies very near 
their initial formation redshift. We can now use the optical 
properties of the radio galaxies to eliminate the low redshift 
component of the distribution. 

Figure |3l shows a Monte-Carlo realization of the dN/dz 
distribution in Figure ^ The shaded histogram represents 
the galaxies that pass our selection criteria, assuming that 
all galaxies have the colors of the solar metallicity instan- 
taneous burst with Zform = 5. These histograms show 
that using optical properties is an effective way of screening 
out the low redshift component of the radio galaxy popula- 
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Fig. 2. — Color-redshift diagrams for two PEGASE models at three formation redsliifts. Tlie upper tracks represent a solar metallicity instantaneous burst model 
with no obscuration, while the lower tracks show a zero initial metallicity exponentially declining star formation with t = 0.5 Gyr and "spheroid" extinction 
assumed. The diagonal lines represent 5 u limits in SDSS g, r, and i bands. Objects above the diagonal hne will not be detected in the individual SDSS band. 



. 1 

- r 

: J 


1 1 1 1 1 1 1 1 _ 






- - 



2 4 6 

KedshifL 



Fig. 3. — Monte-Carlo realization of the dN/dz distribution of FigurefTI 
The shaded histogram represents objects that pass our optical color cuts, as- 
suming that all galaxies have a formation redshift Zform = 5.0. 



tion. The PEGASE models predict that the radio galaxies get 
much bluer and brighter in the optical near their formation 
epoch due to the initial burst of star formation in the mod- 
els. As the histogram in Figure |3] shows, some fraction of 
these galaxies will be excluded from optically selected sam- 
ples near their formation epoch as they become bright and 
blue enough to be detected in most optical surveys. These 
model colors are especially sensitive to the assumptions made 
for the star formation (e.g. obscuration, e-folding time), and 
it is difficult to infer how accurate these color tracks will 
be near the initial starburst. We also note that there is evi- 
dence fo r mode rate obscuration in some high redshift galax- 
ies (e.g..|Dev. So inrad. & Dickinson 1995; Ouchi et al. 200^ 
IVillar-Martin et al...2005i: .Charv et aLi2005) . Nevertheless, 
our models show that we may be somewhat biased against se- 
lecting unobscured galaxies with ongoing star formation near 
their formation epoch. Lowering the optical magnitude cutoff 
would allow these galaxies to enter our sample, although at 
the expense of the low redshift cutoff. Since we are primar- 



ily interested in galaxies with evolved stellar populations even 
at high redshift, as well as maintaining the efficiency of our 
search, we retain the optical selection to favor the low redshift 
cutoff at this minor expense of completeness. 

2.1. Selection using SDSS ant/ FIRST 

We begin with the 2003 April 1 1 version of the FIRST cat- 
alog and select all sources with integrated flux 5i.4ghz > 
100 mJy. The primary motivation for this cut is to de- 
fine a manageable sample size. Flux limits as low as -^lO 
mJy are reasonable for selecting high-redshift radio galaxies, 
and lower flux density limits will be explored in subsequent 
work. The S'i.4ghz > 100 mJy objects are positionally cross 
matched with photometric data from SDSS Data Release 2 
(DR2), which covers more than 3300 square degrees. Ob- 
jects with candidate identifications in any of the u, g, r, ior z 
bands within a conservative radius of 8" are excluded. Radio 
sources with no cataloged SDSS counterpart are visually in- 
spected to further exclude possible low signal-to-noise optical 
counterparts, as well as identifying and excluding extended or 
multiple-lobe radio sources with an obvious optical counter- 
part located some distance from the cataloged radio position 
(an optical source lying between independently cataloged ra- 
dio lobes, for example). Note that we are not restricted to 
unresolved radio sources: Our selection criteria allow for ex- 
tended and multicomponent radio sources to be included as 
well. To further eliminate likely low redshift sources, we op- 
timally combine (Szalav et al. 1999) the g, r, and i SDSS im- 
ages (the three most sensitive of the five SDSS filters) and 
eliminate any radio source with a candidate counterpart in the 
combined image. 

The 5 a point source limiting AB magnitudes for SDSS 
are u = 22 3, g = 23.3, r = 23.1, i = 22.3, z = 20.8, 
ilvezic et al. ' 2000). The coadded g, r, and i SDSS im- 
ages allow us to extend the low redshift range being ex- 
cluded by, in essence, improving our magnitude threshold by 
« 1.0 magnitude. This can be roughly thought of as pro- 
viding an effective 2 a r-band magnitude limit of r « 24.1, 
although the specific limiting value is dependent on the de- 
tails of individual target galaxy SEDs. Coadded images are 
proce ssed with SExtractor version 2.3.2 (Bertin & Arnout^ 
and any objects with a detection greater than 2 a above 
the background are excluded from the sample. Remaining tar- 
get candidates are again visually inspected to eliminate new 
(faint) optical identifications of extended or multi-lobe radio 
sources. This visual inspection introduces some subjectivity 
into the target selection criteria, but it is necessary in order 
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Fig. 4. — Optical color cuts compared to model r — K color as a function 
of redshift. Objects above and to the right of the thick line pass the selection. 



to eliminate obvious SDSS counterparts to double lobed ra- 
dio sources, and sources that are extended in the FIRST cata- 
log. Candidate targets with nearby bright stars or other nearby 
bright confusing sources are also excluded from the final tar- 
get list. These targets are then checked against the NASA 
Extragalactic Database {NEDY to screen out objects that had 
previously been observ ed. Only two objects, 4C —00.62 
(Ikottgering et all Il997h and 3C 257, ( Hewitt & Burbridge 
fl99k .Van Bruegel et alJ fTool were previously identified 
with confirmed redshifts of z = 2.53 and z = 2.474 re- 
spectively. The high redshift of these two radio galaxies sup- 
port the effectiveness of our selection criteria. A final cut of 
RA > 12hr is made to eliminate targets not visible during 
our Spring observing runs. This process yields a total of 172 
target objects. With 2085 S'i.4ghz > 100 mJy radio sources 
in this area of DR2, less than one in ten meet our selection 
criteria, giving us an order of magnitude improvement in effi- 
ciency of finding high redshift radio galaxies over blind spec- 
troscopic targetting of all radio galaxies. 

Figurel^shows a simulation of our multicolor optical selec- 
tion in terms of r — K color as a function of redshift for the 
same two sets of PEGASE spectral synthesis models used in 
FigurelJl The thick line shows the approximate effect of our 
optical selection criteria based on the gri coaddition. Using 
the model dN/ dz of Dunlop and Peacock for radio sources 
brighter than 100 mJy with assumed pure luminosity evolu- 
tion, we expect to exclude all sources with z < 1.7 and an 
average redshift for the sample of z ~ 2.5. 

2.2. Selection on previous surveys 

As a test of our selection criteria, we apply our se- 
lection to subsets of the galaxies in a blind radio survey 

' The NASA/IPAC Extragalactic Database (NED) is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 



jBest et alJ 1 1999ft and the ultr a-steep spectrum (USS) sam- 
ple of IdeBreucket alJ J2001I) . Of the 178 r adio sources 
with spectroscopic redshifts inlBest et al.l Jl999h . only 20 fall 
within SDSS Data Release 3 (DR3). These sources have red- 
shifts ranging from z = 0.004 to z = 2.474, a mean redshift 
of z = 0.78, and a median of z = 0.66. Applying our selec- 
tion criteria leaves only two targets at redshifts of z = 1.339 
and z = 2.474, the only galaxy wit h z > 1.5 in the sa mple. 
Only 9 of the 62 USS sources of de Bre uck et"ai] ("2001) are 
within SDSS DR3, eight of which have firm redshifts. The 
mean and median redshift of this sample are 2.5 and 2.14 re- 
spectively. Applying our selection criteria eliminates the two 
z < 1 galaxies from the sample, as well as one unusually op- 
tically bright galaxy with z = 2.48, where Lya falls in the 
g band, and de Breuck et al. (2001) note that it has relatively 
strong continuum emission. The five remaining galaxies have 
a mean and median of z = 2.58. As these two datasets show, 
our selection criteria is very effective at eliminating the low 
redshift foreground from both blind and USS samples. The 
exclusion of one z > 2 galaxy illustrates the fact that our 
selection may eliminate some galaxies with strong observed 
frame optical emission, and we will explore the extent of this 
bias in a future paper. 

3. THE OBSERVATIONS 

We aquired i^-band images for 96 of the 172 unique can- 
didates on the KPNO 4-meter Mayall telescope over two sets 
of four nights: 2003 April 20-23 and 2004 May 31 to 2004 
June 3, using the Florida Multi-object Imaging Near-IR grism 
Observational Spectrometer (FLAMINGOS) instrument. The 
detector is a 2048x2048 HgCdTe wide-field IR imager and 
multi-slit spectrometer with a pixel size of 0.3165", which 
gives a 10.8' x 10.8' FOV on the 4 meter telescope. Condi- 
tions were highly variable during the 2003 run, with a night 
and a half lost due a combination of high wind and mois- 
ture (first night seeing 1.0" to 1.5", subsequent nights 1.0" 
to 2.3" due to combination of cloud, wind, and moisture). 
The 2004 run was photometric on all four nights with see- 
ing varying from 0.7" to 1.2". Exposures were 20 seconds 
each in the 2003 run and 15 seconds each in 2004 to account 
for brighter sky levels. We observed in a fixed five point 
dither pattern with a separation of 30 arcseconds. Targets 
were observed for 15 minutes unless otherwise noted. The 
5 (7 limiting magnitude in individual images spans K ^ 19.0. 
Three objects among those not detected after fifteen min- 
utes were observed a second time. Objects J1411+0124, 
J1350+0352, and J2242-0808 have total integration times of 
40, 61, and 29 minutes, respectively. Jl 1 23+0530 is the well 
know n z = 2.474 rad io galaxy 3C 257 ( Hewitt & Burbridg^ 
1991; Van Brue gel et al. 1998). one of the most luminous ra- 
dio galaxies known. As it passed our selection criteria, it was 
observed in order to provide a consistency check with previ- 
ous radio galaxy searches. 

3.1. Reduction and Astrometry 

We processed the data using standard NOAO IRAF^ rou- 
tines. A set of dark frames were taken and subtracted from 
each image. A set of 6-10 adjacent (in time) images were 
combined to create a sky flat for subtraction (i.e. "running 

^ IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 
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sky flats"). The images were approximately aligned based on 
the dither offsets, then the IRAF tasks mscgetcat and mscc- 
match were used to accurately register the images for coaddi- 
tion. The I RAF task msccmatch uses a catalog of USN0-A2 
iMonet e t al. 1998) stellar positions and magnitudes for im- 
age registration and transformation, which may include image 
shift, scale change, and axis rotation. The resulting astrome- 
try displayed a systematic offset on most images of between 
0.5" to 1.0", due to sHght differences between the USNO-A2 
and SDSS astrometry. We manually corrected for these sys- 
tematic offsets when constructing catalogs of each field. Our 
final astrometry is accurate to subarcsecond precision, with 
the residual difference between the X-band and SDSS posi- 
tions well fit by gaussian of width 0.25". 

3.2. Photometry 

Each FLAMINGOS field was processed with SExtractor 
version 2.3.2. Quoted i^T-band magnitudes are SExtractor 
MAG_AUTO unless otherwise noted. Due to the large field of 
view of the FLAMINGOS instrument, a large number (■^100) 
of bright sources detected in the 2 Micron All Sky Sur- 
vey ilMASSf' are present in each field. Sources were cross 
matched with point sources from the 2MASS catalog hav- 
ing Kg < 15.3, the 5 a limiting magnitude for point sources. 
Saturated stars were excluded from the comparison. A linear 
least squares fit between the SExtractor and 2MASS objects 
was perfomed to determine the zero point offset between the 
two datasets. This solution was assumed to continue linearly 
beyond the Kg — 15.3 limit. Although we observed in K- 
band, while 2MASS uses Kg, no significant color terms were 
evident (the magnitude comparison between SExtractor and 
2MASS was linear with a slope of unity). For radio galax- 
ies not detected in iiT-band at the 2 a level, we report the 2 a 
magnitude of a point source as a lower limit. 

3.3. Star-Galaxy Separation and ERO definition 

The final ii'-band catalog was positionally cross matched 
with the SDSS catalog using a webservice interface to Open- 
SkyQuery (http : / / www . openskyquery . net) in order 
to obtain r — K colors. A simple nearest neighbor criterion 
was used, and objects with no SDSS counterpart within a 
three arcsecond radius were assigned a 5 cr limiting magni- 
tude ofr = 23.1. 

Star-Galaxy separation was done in two steps: For ob- 
jects detected in SDSS with r < 21.0 we used the SDSS 
star-galaxy classification (probPSF, which is described in 
ISpranton et al. 2002). For fainter r-band objects the separa- 
tion was done by examining the difference between the K- 
band MAG_AUTO magnitude (from SExtractor) and a fixed 
3.5" aperture magnitude returned by SExtractor as a function 
of MAG_AUTO to separate pointlike objects from extended 
galaxies. 

Figure |5l shows the resulting star-galaxy separation. This 
separation was chosen to be in rough agreement with the star- 
galaxy classification from SDSS where reliable classification 
was available. The slope of the resultant number counts for 
bo th stars and galaxi es are in good agreement with those 
of iDaddi et alJ (120001) . Reliable separation was possible to 

' This publication makes use of data products from ttie Two Micron All 
Sky Survey, which is a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/California Institute of Technol- 
ogy, funded by the National Aeronautics and Space Administration and the 
National Science Foundation. 
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green points are r < 21.0 objects classified by SDSS as stars and galaxies 
respectively. Black points above and to the right of the black line are consid- 
ered galaxies. 



K = 17.0, and objects fainter than this are considered to be 
galaxies. 

The exact definition of EROs in the literature is not stan- 
dardized. Some authors use total magnitudes (Wold et alJ 
l2003HCimat ti et al. 2002a), some isophotal/aperture corrected 
JDaddiet alJl200 0). some use matched aperture magnitudes 
in the R- and A'-bands(e.g. Elston, Gonzalez et al. 2006), 
among others. There is the additional variation in the use 
of specific R and K (and Kg) filters used to select EROs, 
making comparison between ERO samples problematic. We 
chose our ERO definition by comparing to publicly avail- 
able data in the Bootes field of the FLAMINGOS Extra- 
galactic (FLAMEX) Survey (Elst on. Gonzalez et al. 2006), 
which lies in the NOAO Deep Wide Field Survey (NDWES)"^ 
( Jannuzi & Dev 1999), as well as SDSS. The Bootes field of 
the FLAMEX survey covers approximately 4.7 square de- 
grees in both J- and i^s-bands. The FLAMEX survey was 
conducted with the FLAMINGOS instrument, and the its lo- 
cation within SDSS allows for a direct comparison with our 
data, although they use the Kg filter where we use K. The 
FLAMEX survey defines EROs based on 6" matched aperture 
magnitudes in NDWFS R and FLAMEX Kg bands. Unfor- 
tunately, the SDSS catalog does not include aperture magni- 
tudes, so we use SDSS model magnitudes in conjunction with 
4" aperture magnitudes in iC-band to define our r — K color, 
as this aperture appears to best mimic the FLAMEX results. 
Cross-matching the FLAMEX, SDSS, and NDWFS data in a 
subset of the Bootes field, we found that an ERO definition of 
^APWB f^'M > 5.50 approxiniated th e R — Ks > 5.0 



definition of lElston Gonzalez et al.l (120061) 



The NOAO Deep Wide-Field Survey (Jannuzi and Dey 1999) is sup- 
ported by the National Optical Astronomy Observatory (NOAO). NOAO is 
operated by AURA, Inc., under a cooperative agreement with the National 
Science Foundation. 
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4. RESULTS 

4.1. The Radio Galaxies 

The results of our X-band observations are shown in Ta- 
bles[2and|2] Tablefjlists objects detected, while Table|2]lists 
non-detections in the infrared. The column entries are: 

1. Object Name: Name of source in lAU J2000 format. 
Objects with double lobe and multiple radio morpholo- 
gies have the infrared object reported first and the indi- 
vidual radio components named a, b, c as necessary. 

2. Date: Date observed. 

3. Radio RA and DEC: Position of radio detection from 
the FIRST catalog (J2000) 

4. iC-band RA and DEC: Position of the infrared detected 
object (J2000) 

5. Extended: Y if object shows visible extension in the 
FIRST image or if the object is a double or multiple 
lobe radio source. N if object is unresolved in the radio 
image. 

6. /peak: Peak 1.4 GHz flux density from FIRST catalog 
(mJy) 

7. fiat- Interated 1.4 GHz flux density from FIRST cata- 
log (mJy) 

8. K: apparent magnitude and error for infrared object. 
Magnitudes are in the 2MASS system, which is slightly 
offset from the CIT system, given by: ii'5(2MASS) ^ 
K(CIT) + 0.0000 X (J - K)ciT - 0.024 ± 0.003 
(Carpenter 2001). For nondetections (Table |3 the 2 a 
detection limit is reported. Errors are those output by 
SExtractor. 

9. Seeing: The full width half maximum seeing measured 
for the field, given in arcseconds. 

10. zk~z' estimated redshift based on a fit to the Hubble 
K — z relation, given hy K — 4.62 log(2;) + 17.2. For 
nondetections (Table|2ji, the estimated redshift is given 
as a lower limit for the 2 a limiting magnitude. 

We observe a total of 96 unique objects, detecting 70 of 
these at a 2a level to a magnitude limit of K « 19.5 in 
most cases. Three objects, J1102+0250, J1237+0135, and 
J1606+4751 show highly extended double radio lobes with 
complex morphologies. Subsequent visual inspection of the 
coadded gri SDSS images for these objects reveals a probable 
optical counterpart that was missed in the original target selec- 
tion. They thus do not meet our selection criteria, and we in- 
clude them here merely for completeness. The total areal cov- 
erage of the A'-band fields (exluding the field of J0742+3256, 
where a transient problem with the detector corrupted a por- 
tion of the field, but did not affect the radio galaxy target) 
is 8570 square arcminutes (2.38 square degrees). Of the 96 
target objects, 27 are double lobed or multiple radio sources, 
having multiple components in the FIRST catalog. All but 
four of these multiple sources are identified in the NIR. Ten 
additional radio galaxies appear to be resolved in the radio 
postage images, and eight of these are identified with A'-band 
counterparts. 



Figure |6l shows images for a selection of our K- 
band detected objects. FIRST contours are overlaid 
with the outer contour at 5mJy, and each subsequent 
contour indicating an increase by a factor of two. 
Postage stamp images for all galaxies can be found at 
http : //lahmu .phyast .pitt . edu/dragons/. 

Figure shows the distribution of if -band magnitudes for 
our observed objects. The solid histogram represents 2 a 
lower limits for our non-detections. The solid vertical line 
represents the magnitude of an object at our expected red- 
shift cutoff of z = 1.7, given our fit to the K — z relation. 
The shaded region represents the approximate range of mag- 
nitudes for objects that follow the K — z relation at this red- 
shift. 

Figure |8l shows the estimated redshift distribution based on 
our linear fit to the K — z Hubble diagram. Caution is advised 
in interpreting this histogram, because of the large intrinsic 
scatter in the K — z diagram. The shaded portion of the his- 
togram represents the lower limits on the estimated redshifts 
for the non-detected objects. The Monte-Carlo simulation of 
the Dunlop and Peacock dN/dz convolved with our selec- 
tion criteria from Figure |3l is shown for comparison. If we 
assign redshifts to the radio galaxies based on the linear fit 
to the K — z diagram, the mean redshift for this sample is 
z = 2.5 and the median redshift is z = 2.0. Comparison 
with the expected redshift distribution in Figure |S] shows that 
we have 35 objects with K < 18.3 that we expected to be 
excluded from our sample based on our model color tracks. 
While three of these objects are resolved in the radio images 
and subsequent closer examination revealed counterpart ob- 
jects in SDSS, as indicated above, a full third of our sam- 
ple are brighter in K than expected. This can be partially 
attributed to the scatter in the K — z relation if these galaxies 
are at the upper end of the redshift range given their K mag- 
nitude. For example, the previously observed radio galaxy 
3C 257 (J1123H-0530) is known to be at z = 2.474, while 
its K — 17.51 magnitude corresponds to z ~ 1.2 on a lin- 
ear fit to the K — z diagram, indi cative of the large intrin- 
sic scatter in the K — z relation. IVan Bruegel et alJ J1998I) 
point out that Ha falls into the K-hand for 3C 257, which 
may explain its relative brightness in K. Similar line con- 
tamination may be responsible for some of the brighter than 
expected objects in our sample, but even this does not fully 
explain the brightest and reddest of our objects. Thirteen ob- 
jects with K < 17.5 have anomalously red g — K,r — K, and 
i — K colors that are not fit by even our extreme model tem- 
plates. Galactic r-band extinction from the SDSS database 
for these objects are listed in Table |3] based on the dust 
maps of ISchlegel. Finkbeiner. & DavisTj 19981) . There is a 
modest amount of extinction for three objects: J0941H-0127, 
J1548H-0036, and J1604-0013, though not enough to fully ex- 
plain their extreme color. Applying the extinction model of 
Calzetti et al. (2000) to the non-evo lving elliptical template 
of Colema n. Wu. & WeedmanI ( fT98QD and assuming a red- 
shift of near unity given by the K — z diagram, we require 
extinctions of Ay > 0.5 — 1.5 in order to reach the lower 
limit r — K colors of these infrared bright sources. Determin- 
ing spectroscopic redshifts for these objects will be the best 
way to determine whether the anomalously red colors are due 
to obscuration at lower redshift, a more luminous galaxy at 
higher redshift, or some combination of the two. These ob- 
jects could also be a separate class not covered by our galaxy 
models, e.g. high redshift radio loud quasars or low redshift 
type II quasars where the AGN is completely obscured. 
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Fig. 6. — A selection of if-band images. North is to the right and East is down. FIRST radio contours are overlaid with an outer contour at 5mJy, and each 
subsequent contour indicating an increase by a factor of two. 





2 4 
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Fig. 7. — Distribution of ft'-band magnitudes for the observed objects. The 
filled histogram represents the 2 cr limiting magnitudes of the non-detections. 
The solid vertical line represents the A'-band magnitude for a galaxy at our 
expected lower redshift limit of z = 1.7 and the shaded area represents the 
approximate scatter given the K — z diagram at this redshift. 



Fig. 8. — Estimated redshift distribution based on a linear fit to the K — z 
Hubble diagram. The shaded histogram represents lower hmit redshift es- 
timates from the magnitude lower limits. The dashed histogram sh ows the 
distribution expected from convolving the lPunlop & PeacocU1199(D dN/dz 
with our selection criteria. 
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Fig. 9. — Spectral slope as a function of zk-z for the 75 targets with 
Saes MHz flux densities measured in the Texas survey. The horizontal line 
represents the USS cut of Van Bruesel et al. 1 1998) 



4.2. Radio Spectral Index 

Table@] shows the radio spectral index (a, given by Sy oc 
v^") computed from FIRST and the Texas 365 MHz survey 
((bouglas et al. 1996). Van Bruegel et al. ( 1998) used an ultra- 
steep spectrum (USS) cut of > 0.8 to select a sample 
of high redshift radio galaxy candidates. Of our 96 targets, 75 
have observations at 365 MHz, and of these 24 have a < 0.8 
and woul d not be selected 

by the iVan Bniegel gtlTI llT998l) criteria. Figure |9l shows 
the spectral 

index of our targets as a function of zk-z- FigurelTolshows 
a histogram of the redshifts inferred from the K — z diagram 
for two spectral index ranges, flatter and steeper than a — 
0.8. The median redshift of the a < 0.8 sample is zk-z — 
2.0 and the mean is zk-z — 2.13, compared to a median of 
zk-z ~ 2.0 and mean zk^z = 2.17 for the total sample, 
and median zk-z ~ 2.0 and mean zk-z — 2.30 for the 
a > 0.8 sample. This suggests that flatter spectral slope 
systems may contribute a non-negligible fraction of the HzRG 
population. The insensitivity of our method to spectral slope 
allows us to select candidates that would be missed by USS 
selection techniques. Applying an USS criteria to our dataset 
would eliminate a full third of the targets. Figure[^shows the 
1.4 GHz luminosities of the 75 galaxies given their aggg^ and 
K — z redshift. The three lines represent 100 mJy, 500 mJy, 
and 1 Jyl.4 GHz flux limits, assuming a = 0.4. 

Table^Jlists a subset of our candidates with additional radio 
observations at 4.85 GHz (From Gregorv & Condon I99w 
and 151 MHz (From the 6C survey. Hales et al. 1988, 199^- 
This table shows the frequency dependence of the two point 
spectral slope. Several of the sources show significant devia- 
tions from a power law over the frequency range in question, 
showing that the USS sample selection will differ depending 
on frequencies used, while our selection method is unaffected 
by objects with concave radio spectra. 
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Fig. 10. — Histogram of redshifts infen'ed from the K — z diagram divided 
into a < 0.8(shaded) and ultra-steep a > 0.8 emphasising the presence at 
high redshift (z/f — z > 2) of flatter spectrum sources. 
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Fig. 1 1 . — 1.4 GHz radio luminosity as a function of K — z redshift. The 
lines represent 100 mJy (solid), 500 mJy (dashed), and 1 Jy (dot-dashed) flux 
density limits assuming a spectral slope a = 0.4. Squares mai'k objects with 
OggI'' > 0.8, triangles indicate Qagg < 0.8. The solid black circle marks 
the radio luminosity of 3C 257 witli spectroscopic redshift z = 2.474. 



4.3. Environment 

If, as we assert, high redshift radio galaxies form in the 
most overdense regions of the early Universe, then hierar- 
chical formation scenarios predict an enhanced number of 
(proto-)galaxies associated with this overdensity. Although 
the strong spatial clustering of EROs is often attributed 
to association with high redshift galaxy overdensities, it is 
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Fig. 12. — Diffe rential A"-band galaxy number counts . For comparison the 
numbe r co unts of [ G ardner etal] 11993l) : ISzokolv et alj 11993) : llbtani etalj 
O051 and lCimatti et alJI2002U) ai'e shown. 



only recently that direct evidence for this has been found 
dOeorgakakis et al. 2005). We search for overdensities of both 
EROs and isT-band selected galaxies in the vicinity of our 
radio galaxy candidates. We begin with a consideration of 
the source counts of both X-band objects, as well as EROs. 
Figure [21 shows the A'-band differential number counts for 
the galaxies and stars in our 95 fields. No completeness 
corrections have been app lied. Plotted for comparison are 
counts from lGardner et a llJl993h.lSzokolv et al.l(fl998) (A%), 
iTotani et al. ( 2001), and Cimatti et alJ (I2002at~ rom the A:20 
survey, also Kg). The raw counts are given in Table|6l As 
the Figure shows, we are in good agreement with previous in- 
frared observations, and it appears that the radio galaxy target 
selection does not bias our AT -band counts to the magnitude 
limits probed. 

Figure ^] shows the r — AT as a function of K color- 
magnitude diagram after cross-matching our infrared catalog 
with the SDSS. We define an Extremely Red Object (ERO) 
as having r — AT > 5.50 as indicated by the horizontal line 
on this Figure. Objects having no corresponding SDSS object 
within 3" are assigned a limiting magnitude of r = 23.1. This 
yields 

a total of 479 EROs, and Figure shows the cumulative 
ERO surface density for this sample. Shown for comparison 
are the SDSS/FLAMEX ERO cumulative density for the « 
4.7 degrees^ of the Bootes field. As can be seen, our ERO sur- 
face density for AT > 17.0 EROs is higher than the surface 
densities of both Daddi et al. (2000) and the FLAMEX/SDSS 
data, suggestive of an ov erdensity of EROs associated with 
our radio galaxy targets. iWold et aP (j2003) observe EROs 
surrounding z ~ 2 radio loud quasars, and the Figure shows 
a higher ERO surf a ce density than the "fiel d" surveys of 
[Roche et al. (2002), Dad di et alJ |2000), and 'Cimat ti et aP 
(|2002a). Our data seem to show a continuation of the trend 
seen by Wold et al. (2003). Because of the r = 23.1 limit 
from SDSS we are limited to studying AT < 17.60 EROs. Due 



Fig. 13. — r — if as a function of K, showing star-galaxy separation 
and ERO definition. The diagonal upper limit in the measurements is due 
to the 5 cr r-band limit from the SDSS assigned to non-detections, while the 
horizontal line marks the r — K = 5.50 color used to define EROs. Red 
points are objects classified as stars, while black are galaxies. The diagonal 
line represents r = 21.0, below which the SDSS star-galaxy separation is 
used. No star-galaxy separation is done for objects above this diagonal line 
and to the right of the vertical line at iC = 17.0. Open circles represent 2 a 
r — K lower limits for the target radio galaxies in these fields. 
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Fig. 14. — Cumulative ERO surface density for our data set and 
FLAMEX/SDSS. Sh o wn for comparison are data fromlDaddi et alj t200a) : 
ICimatti et all 12002 J) : lRoche et alJ<20()3) and lWold et alJ 120031) . 



to the small areas covered by previous ERO studies, surface 
densities at this bright K magnitude are uncertain. 5.9% (30 
out of 509) o f our EROs are classifi ed a s stars, smal l er than 
the « 9% of lManucci et'aP ll2()0l and iWold et alJ ll2()0l . 
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Fig. 15. — Average surface density of K < 17.50 EROs as a function of 
radial distance from the radio galaxy for 95 target fields. The four panels are 
subsamples selected by the A'-band magnitude of the target radio galaxy. The 
horizontal line represents the "local field density" of the sample, defined as 
the average surface density of EROs between 120" and 180" from the radio 
source. 
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Fig. 16. — Surface density of K < 19.0 galaxies as a function of radial 
distance from the radio galaxy. The four panels are subsamples selected by 
the _R'-band magnitude of the target radio galaxy. The horizontal line in each 
panel represents the "local field density" of _fC-band selected galaxies, defined 
as the average surface density of galaxies between 120" and 180" from the 
radio source. 



which may be a concern. At these bright if -band magnitudes, 
robust star galaxy separation is increasingly important. We 
note, again, that our overall star and galaxy counts are in good 
agreement with Dad di et all (|2000) in the relevant magnitude 
range, leading us to believe that our star-galaxy separation 
is robust. Figure \^ shows that our A'-band number counts 
are consistent with previous surveys, yet our ERO density is 
higher at K > 17.0 than in the same surveys. This is because 
only a small percentage of the X-band selected galaxies are 
EROs, and the increase in this rare population is small com- 
pared to the overall galaxy counts. Thus, the EROs appear to 
trace an overdensity not readily apparent in the A'-band data 
alone, as pointed out by Wold et al. (2003). 

We now examine the clustering of EROs around our radio 
galaxies. We restrict our analysis to K < 17.50 EROs in 
order to remain complete given the r\nn = 23.1 magnitude 
limit of our SDSS data. Because we expect the radio galaxy 
to be, by far, the most luminous galaxy in its local environ- 
ment, we expect to see a possible overdensity around only the 
brightest of our targets. An overdensity of EROs with A'-band 
magnitudes brighter than the radio galaxy would be a possi- 
ble indication of a foreground structure. Field by field, the 
density of EROs is quite inhomogeneous, as expected. More 
than half of the fields have no bright EROs within 100" of 
the radio galaxy target, while seven fields have three or more 
EROs within this radius. The seven densest fields have ra- 
dio galaxy targets with magnitudes evenly spaced in the range 
K = 17.4 — 20.6, showing no trend of high ERO density 
with radio galaxy NIR brightness, contrary to our expecta- 
tions. Figure^lshows the average distribution of EROs with 
K < 17.50 as a function of radial distance from the radio 
galaxy for 95 of our target fields. We do not include the ra- 
dio galaxy itself, although some are EROs. The horizontal 
line is a measure of the "local field" density of EROs, defined 



as the density of EROs between 120" and 180" from the ra- 
dio source. There is a modest overdensity of EROs within an 
arcminute of the radio targets , although the uncertainties are 
very large. IWold et all J2003') perform a similar analysis and 
see no such evidence for clustering for EROs brighter than 
Ks — 19.5 around their z ^ 2 radio loud quasars. A deeper 
sample of EROs is necessary for a complete analysis of the 
radio galaxy environment. This could be established from our 
existing data through the addition of deep i?-band imaging 
(R « 25), and we are actively pursuing this goal. Deep R- 
band imaging would also allow for more rigorous star-galaxy 
separation to fainter magnitudes, to better determine the stel- 
lar contamination of our ERO population. 

Figure [T6l shows the distribution of AT < 19.0 galaxies 
around our targets. The four panels of Figure[^represent the 
division of the sample by the apparent magnitude of the target 
radio galaxy. Given the K ~ z Hubble diagram, this is a rough 
proxy for redshift. The horizontal line in each panel shows a 
measure of the "local field" density of galaxies (defined as 
the density of objects between 120" and 180" from the radio 
source). No clustering is seen in any of the bins. This is not 
unexpected, as .Hall & Green (1998) and Hall et al. ( 2001 ) do 
not see an excess around z ^ 1.5 radio loud quasars below 
K — 19, but do see an excess of galaxies at K = 19 — 
20.5. Deeper AT-band imaging is necessary to detect possible 
companion galaxies at the high redshifts of our radio galaxy 
targets. Such imaging will be conducted on future additions 
to our target sample (see §|5lfor details). 

5. SUMMARY AND FUTURE WORK 

Using a novel radio-optical selection technique, we have 
obtained NIR imaging for 96 HzRG candidates. Based on red- 
shift estimates from the K—z Hubble diagram, our new selec- 
tion technique appears to be effective at identifying HzRGs, 
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as more than half of our targets have if -band photometry con- 
sistent with z > 2. This technique is not sensitive to ra- 
dio spectral slope, and avoids the frequency dependence of 
USS techniques for galaxies with non-power law radio spec- 
tra. Of the 75 target galaxies present in the Texas 365 MHz 
survey 24 have al'g^Rm ^ < 0.8. These would be excluded 
by the USS criteria of VanB rueeel et alJ jl998h . and all but 
four of these 75 have asgs'j^^ < 1-3, whi ch would be ex- 
cluded from more recent USS selections (e.g.'de Breuck et al] 
I2OOI: Jarvis et al. 2001). Comparison with the previous sur- 
veys of Best et al. ( 1999) and de Breuck et al. (2001) shows 
that our technique selects nearly all of their high redshift tar- 
gets, while also eliminating the low redshift sources. The K- 
band number counts are consistent with previous work, and 
the large non-contiguous area (2.38 square degrees) covered 
by the DRaGONS survey to date makes it an excellent re- 
source for exploring galaxy properties through the combina- 
tion of NIR and optical data. 

We have uncovered a previously unseen class of radio 
sources with anomalously red colors (r — K > 6.5 — 7), 
which may be evidence of significant obscuration at moderate 
redshifts. These galaxies represent more than ten percent of 
our observed sample, indicating that they may be a substan- 
tial percentage of the radio galaxy population. Being non- 
detections in SDSS (with r > 24.1), cuiTent optical AGN 
selection techniques are insensitive to these sources. Radio 
loud QSOs comprise only ~ 5 — 10% of overall QSOs; there- 
fore, these objects could represent a significant contribution 
to radio loud AGNs that are not counted in current samples. 

We have identified a sample of 479 EROs to a depth of K = 
17.50, one of the largest bright ERO samples to date, and have 
identified a modest overdensity of EROs around our HzRG 
candidates. We do not detect any similar overdensity when all 
the if -band selected galaxies to K — 19.0 are considered. 

Spectroscopic follow up of the radio galaxy candidates and 
the surrounding EROs is of the utmost importance, both for 
validation of our HzRG candidates, as well as determination 
of the relationship between the radio galaxies and the ERO 
overdensity. Firm redshifts will also allow us to accurately 
test the z — a relation for our non-USS selected sample. Near 
infrared spectroscopy of the anomalously red radio galaxies 
should pinpoint the cause of their extreme colors, as well as 
determine the amount of obscuration present at intermediate 
and high redshifts. 

Deeper optical imaging is required in order to establish col- 
ors for the radio galaxies and explore their morphology in 
more detail. Deeper optical imaging will also allow us to se- 
lect a deeper and more complete sample of EROs, allowing 
for a more detailed study of the radio galaxy environment and 
the nature of the ERO overdensity. A more quantitative study 



of ERO clustering will also be possible with deeper optical 
data (including a two-point angular correlation function anal- 
ysis). 

A future paper will also examine the effects of both our ra- 
dio flux density cut, as well as our optical color cuts. In future 
observing runs we plan to target a small sample of lower flux 
density radio sources that meet our optical color cuts. We will 
supplement this sample with galaxies from the FLAMEX sur- 
vey that meet our selection criteria. This dataset will allow us 
to examine any biases introduced by our radio selection. Sev- 
eral of the targets in our most recent round of observing lie in 
the SDSS Southern Survey, an area of SDSS that is more than 
four times more sensitive than the standard SDSS survey. This 
deeper imaging data will allow us to study the optical proper- 
ties of those targets that are detected. 
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-06:03:00.34 


N 


20:59:32.87 


-06:02:59.59 


155.63 


161.72 


18.32±0.10 


0.9 


1.7 


J2107^0701 


02 June 2004 


21:07:45.46 


-07:01:07.83 


N 


21:07:45.47 


-07:01:06.54 


523.60 


550.60 


18.37±0.09 


0.9 


1.8 


J2223-0757 


01 June 2004 


22:23:26.52 


-07:57:08.07 


N 


22:23:26.51 


-07:57:06.94 


104.79 


108.34 


17.70±0.08 


0.9 


1.3 


J2247-0910 


01 June 2004 


22:47:23.79 


-09:10:49.74 


Y 


22:47:23.69 


-09:10:49.26 


77.10 


103.37 


17.99±0.09 


0.9 


1.5 



TAB LEI — Continued 



Object Name" 



date 



radio RA 



radio DEC Extended? 



KRA 



KDEC 



/pcak(mJy) /int(mJy) 



K 



seeing( ) 



ZK- 



J2309-0834 
J2316-0846 
J2336-0838 
a 
b 

J2337-0852 



03 June 2004 

03 June 2004 

04 June 2004 



04 June 2004 



23:09:04.29 
23:16:35.08 

23:36:18.11 
23:36:18.65 
23:37:32.44 



-08:34:57.19 
-08:46:17.72 



-08 



38:43.58 
38:48.27 
52:39.49 



N 
N 
Y 



23:09:4.29 
23:16:35.07 
23:36:18.35 



-08:34:56.83 
-08:46:17.16 
-08:38:48.87 



23:37:32.46 -08:52:38.83 



103.19 
101.85 

124.78 
104.16 
118.04 



111.67 
105.40 

134.50 
110.32 
120.14 



18.23±0.12 
18.59±0.14 
18.50±0.15 



18.20±0.13 



0.9 
0.9 
0.9 



1.7 
2.0 
1.9 



1.6 



Objects named with u and b are double lobed sources with the radio properties hsted corresponding to the object above them on the table 
Based on a fit to the - z Hubble diagram it = 4.62 * log(z) + 17.2 

This object is an extended double lobe with a corresponding faint match in SDSS and should not have been included in our sample. The object has SDSS magnitudes g=23.25 r=22.29 i=21.26 z=19.90 

This object, 3C 257 was previously identified as a z — 2.474 radio galaxy ^^^^^^^^^^^^^^J. The KT-band contains Ha emission, leading to the underestimate of redshift from the K-z diagram (set^^^^^^^^^^^^^ for details). 
This object is a wide separation double lobe with several candidate objects along the radio axis. The probable match is faintly delected in SDSS with magnitudes org-22.6 r=22.13 i-21.32 z=20.27 
Object J141 1+0124 was observed on subsequent nights for a total integration time of 40 minutes 

This object has a match in SDSS g-22.86 r=22.45 i=22.11 z-21.12, was improperly targeted because it is a weak double lobe 
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TABLE 2 

Non-Detections 



Object Name 


date 


radio RA 


radio DEC 


Extended? 




/peak(mJy) 


/int(mJy) 


K (2 a) 


seeing( ) 




J1022+0357 


23 April 2003 


10:22:01.03 


+03:57:37.56 


N 


194.31 


200.86 


>20.11 


2.3 


>3.2 


J1028+0144 


21 April 2003 


10:28:02.79 


+01:44:06.51 


Y 


56.82 


114.81 


> 19.80 


2.0 


>2.7 


J 1044+05 3 8 


20 April 2003 


10:44:19.88 


+05:38:07.98 


N 


125.54 


130.20 


>20.32 


1.2 


>3.5 


J 1047+02 16 


23 April 2003 


10:47:11.32 


+02:16:28.21 


N 


122.32 


125.81 


>20.09 


2.2 


>3.1 


Jl 144+0254 


23 April 2003 


11:44:34.26 


+02:54:25.56 


N 


114.01 


119.18 


>20.13 


2.1 


>3.2 


J 122 1+0248 


23 April 2003 


12-21 -39 93 


+02:48:28.01 


N 


139.93 


142.87 


>20.13 


1.9 


>3.2 


J 1234+0024 


03 June 2004 


12:34:30.79 


+00:24:59.45 


N 


94.84 


100.55 


>20.35 


1.2 


>3.6 


J1240— 0017 


21 April 2003 


12:40:12.23 


—00:17:30.34 


N 


137.73 


150.90 


>20.28 


2.0 


>3.4 


J1303+0026 


23 April 2003 


13:03:57.48 


+00:26:45.41 


N 


91.23 


104.04 


>20.15 


2.0 


>3.2 


11314+0330 


21 April 2003 


13:14:22.82 


+03:30:22.14 


Y 


163.59 


250.05 


>20.31 


1.6 


>3.5 


J 1329+0 133 


01 Jvine 2004 


13:29:18.78 


+01:33:40.80 


N 


87.80 


102.64 


>20.38 


1.3 


>3.6 


11350+0352" 


2003-2004 


13:50:24.37 


+03:52:43.90 


N 


99.51 


104.28 


> 21.02 


1.1 


>5.0 


J1421+0248 


01 June 2004 






Y 






>20.33 


1.0 


>3.5 






14:21:10.957 


+02:48:35.76 




152.15 


165.01 








b 




14:21:11.20 


+02:48:29.10 




158.20 


177.62 








J1431+0511 


04 June 2004 


14:31:09.58 


+05:11:17.85 


Y 


1.31 


2.22 


>20.19 


0.9 


>3.3 






14:31:08.10 


+05:11:21.18 




216.42 


226.17 








b 




14:31:10.86 


+05:11:15.80 




158.20 


177.62 








J 1500+0031 


23 April 2003 


15:00:55.34 


+00:31:58.52 


N 


141.31 


145.90 


>20.21 


2.0 


>3.3 


J1507+6003 


04 June 2004 


15:07:44.31 


+60:03:12.68 


N 


181.08 


189.23 


>20.18 


0.9 


>3.3 


J1527+4352 


02 June 2004 


15:27:51.49 


+43:52:4.79 


N 


138.95 


140.11 


>20.35 


0.9 


>3.6 


J1554+3942 


03 June 2004 






Y 






>20.23 


0.9 


>3.3 


a 




15:54:17.02 


+39:42:27.44 




35.95 


55.94 








b 




15:54:17.81 


+39:42:18.95 




98.92 


113.29 








J 1557+4657 


04 June 2004 


15:57:24.61 


+46:57:54.29 


N 


198.10 


204.20 


>20.14 


0.9 


>3.2 


J1618+5210 


04 June 2004 


16:18:55.60 


+52:10:41.40 


N 


108.26 


112.29 


>20.21 


0.9 


>3.3 


J1641+4209 


22 April 2003 


16:41:30.14 


+42:09:25.99 


N 


260.20 


263.68 


>20.37 


1.0 


>3.6 


J1648+3623 


04 June 2004 






Y 






>20.12 


0.9 


>3.2 


a 




16:48:51.58 


+36:23:39.10 




121.33 


128.45 








b 




16:48:53.02 


+36:23:24.56 




123.58 


128.96 








J1700+3830 


20 April 2003 


17:00:19.95 


+38:30:33.93 


N 


428.65 


430.20 


>20.40 


1.2 


>3.6 


J1711+3047 


21 April 2003 


17:11:26.65 


+30:47:45.89 


N 


124.27 


124.77 


>20.41 


1.2 


>3.7 


J2221-0901 


02 June 2004 


22:21:48.04 


-09:01:58.95 


N 


224.64 


234.50 


>20.32 


0.9 


>3.5 


J2242-0808'' 


01 & 04 June 2004 


22:42:34.05 


-08:08:21.86 


N 


125.64 


129.97 


>20.63 


0.9 


>4.1 



J135(>f0352 was observed both 20 April 2003 and 02 June 2004, for a total integration time of 61 minutes with no detection 
J2242 — 0808 was observed on two nights in 2004 for a total integration time of 29 minutes 



TABLE 3 

Extinction FOR K < 17.5 Galaxies 



Name 


K 


Extinction 


0941+0127 


17, 


,17 


0.32 


0958+0324 


17, 


,46 


0.08 


1155+0305 


17, 


,36 


0.10 


1208+4943 


17, 


,03 


0.06 


1332+0101 


17, 


.05 


0.08 


1400+0053 


17, 


.23 


0.10 


1408+0116 


17, 


.08 


0.11 


1423+0139 


17, 


.37 


0.09 


1532+4432 


17, 


.16 


0.05 


1541+5259 


17, 


.24 


0.03 


1548+0036 


16, 


,83 


0.23 


1604-0013 


17, 


,13 


0.34 


1649+3350 


16 


.86 


0.06 



TABLE 4 
Radio Spectral Index 



Name 


FIRST 1.4 GHz (mjy) 


Texas 365 MHz (mJy)' 


„1400 
■^365 


J0742+3256 


144.60 


426.0 


0.80 


J0831+5210 


246.00 


837.0 


0.91 


J0941+0127 


120.56 


290.0 


0.65 


J0958+0324 


638.00 


2666.6 


1.06 


J 1022+0357 


200.86 


823.0 


1.05 


J1028+0144 


114.81 


391.0 


0.91 


J1044+0538 


130.20 






J1047+0216 


125.81 


311.0 


0.67 


Jl 102+0250 


161.00 


1626.0 


1.72 


Jl 123+0530 


1743.00 


5903.0 


0.91 



DRAGONS 



TABLE 4 — Continued 



Name 


FIRST 1.4 Otlz (mjy) 


Texas 365 MHz (mJy) 


1400 
"366 




212.60 


982.0 


1.14 


Jl 144+0254 


119.18 


445.0 


0.98 


11 1 55+0305 


162.80 


591.0 


0.96 


J1208+0414 


641.00 


2212.0 


0.92 


T 1208+4943 


198.98 


678.0 


0.92 


J1221+0248 


142.87 


510.0 


0.95 


11934+0024 


100.55 


700.0 


1.44 


J1236+0150 


154.53 


558.0 


0.96 


J 1237+0 135 


426.00 


1971.0 


1.14 


J 1240-00 17 


150.90 


419.0 


0.76 


11250+6043 


304.05 


725.0 


0.65 


Tl 25Q+055Q 


178.90 


806.0 


1.12 


Tl 303+0026 


104.04 


401.0 


1.00 


Tl 308-0022 


241.17 






11312+0009 


112.60 


490.0 


1.09 


Tl 3 13+6250 


132.18 


476.0 


0.95 


Tl 3 14+0330 


250.05 


845.0 


0.91 


11315+0533 


146.96 


423.0 


0.79 


11329+0133 


102.64 






J1332+0101 


411.00 


1430.0 


0.93 


J 1336+0207 


127.76 


410.0 


0.87 


J1350+0352 


104.28 


546.0 


1.23 


J 1400+005 3 


130.80 






J 1402+0342 


540.40 


1192.0 


0.59 


J 1403+6048 


794.62 


1946.0 


0.67 


J1408+0116 


612.75 


1325.0 


0.57 


J141 1+0124 


186.50 


1029.0 


1.27 


J 142 1+0248 


342.63 


984.0 


0.78 


11423+0139 


211.65 


394.0 


0.46 


J1431+0511 


226.17 


956.0 


1.07 


Tl 438+01 50 


118.50 


413.0 


0.93 


J1438+6149 


121.27 






J1451+5404 


554.63 


2108.0 


0.99 


J1452+0032 


639.00 


3133.0 


1.18 


Tl 500+0031 


145.90 


549.0 


0.99 


Tl 507+6003 


189.23 






J15 10+5244 


505.76 


1596.0 


0.85 


J15 15+5744 


144.16 


400.0 


0.76 


J1523-0018 


231.44 


444.0 


0.48 


J1526+0408 


156.40 


504.0 


0.87 


J1527+4352 


140.11 






J 1532+4432 


249.65 






J 154 1+5259 


193.00 


565.0 


0.80 


J1543+571 1 


103.39 


279.0 


0.74 


J1547+4839 


214.75 






J 1548+0036 


126.00 


313.0 


0.68 


1154S-0033 


432.00 


1221.0 


0.77 


J1549+4719 


106.33 


252.0 


0.64 


Tl 554+4729 

J X^^ 1 1 1 / ^y 


149.73 


612.0 


1.05 


Tl 554+3942 


169.23 


697.0 


1.05 


Tl 557+4657 


204.20 


756.0 


0.97 


Jl 559+501 1 


130.28 


419.0 


0.87 


J1604+4746 


366.10 






J 1604-00 13 


128.28 






J 1606+4751 


109.02 


670.0 


1.35 


J 1609+3700 


102.52 


206.0 


0.52 


J1617+4848 


244.32 






J1618+5210 


112.29 


294.0 


0.72 


J 1629+4937 


197.05 


499.0 


0.69 


J1632+4056 


203.15 






J1634+4155 


263.47 


1155.0 


1.10 


J1636+4808 


263.30 


981.0 


0.98 


J1637+3223 


147.47 


438.0 


0.81 


J1641+4209 


263.68 


1186.0 


1.12 


J1643+4518 


110.39 


304.0 


0.75 


J1645+4152 


115.96 






J1648+4233 


169.30 






J1648+3623 


257.00 


1197.0 


1.14 


J1649+3350 


163.91 






J1654+4125 


228.56 


600.0 


0.72 


J1655+2723 


169.53 






J 1656+2707 


164.76 


347.0 


0.55 


J1700+3830 


430.20 


868.0 


0.52 


J1707+2408 


169.15 


787.0 


1.14 


J171 1+3047 


124.77 


407.0 


0.88 


J1715+3027 


385.34 


1108.0 


0.79 
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TABLE 4 — Continued 



Name FIRST 1.4 GHz (mJy) Texas 365 MHz (mJy)' a. 



J2059-0603 


161.72 






J2107-0701 


550.60 


2373.0 


1.09 


J222 1-0901 


234.50 


930.0 


1.02 


J2223-0757 


108.34 


361.0 


0.90 


J2242-0808 


129.97 


573.0 


1.10 


J2247-0910 


103.37 






J2309-0846 


111.67 






J23 16-0846 


105.40 


338.0 


0.87 


J2336-0838 


244.82 


1137.0 


1.14 


J2337-0852 


120.14 


507.0 


1.07 



1 From lDouglas et alj < 19961) 



TABLE 5 
Radio Spectral Index 



Name 


4.85 GHz (mJy)' 


FIRST 1.4 GHz (mJy) 


Texa.s 365 MHz (mJy) 


151 MHz (mJy) 


1400 
"365 


„4850 
"1400 


1400 
"l51 


365 

"151 


J1123-H0530 


526 ± 73 


1743.0 ± 69.7 


5903 ± 140 




0.9 ± 0.041 


1.0 ±0.12 






J 1208+4943 


58 ± 8 


199.0 ± 8.0 


678 ± 34 


990 ± 103 - 


0.9 ± 0.053 


1.0 ±0.12 


0.7 ± 0.052 


0.4 ±0.13 


J1313+6250 


29 ±6 


132.2 ±5.3 


476 ± 19 


970 ± 48.5^ 


1.0 ±0.048 


1.2 ±0.17 


0.9 ± 0.032 


0.8 ± 0.072 


J 1403+6048 


268 ± 26 


794.6 ± 31.8 


1946 ±31 


3070 ± 154' 


0.7 ± 0.039 


0.9 ± 0.088 


0.6 ± 0.032 


0.5 ± 0.060 


J 1543+57 11 


25 ±5 


103.4 ±4.1 


279 ± 15 


490 ± 24.5^ 


0.7 ± 0.055 


1.1 ±0.17 


0.7 ± 0.032 


0.6 ± 0.083 


J 1554+472 


34 ±6 


149.7 ± 6.0 


612 ± 18 


1260 ± 126- 


1.0 ±0.043 


1.2 ±0.15 


1.0 ±0.050 


0.8 ±0.12 


J1557+4657 


39 ±6 


204.2 ± 8.2 


756 ± 19 


1580 ± 158^ 


1.0 ±0.042 


1.3 ±0.13 


0.9 ± 0.050 


0.8 ±0.12 


J1559+5011 


37 ±6 


130.3 ±5.2 


419 ± 19 


960 ± 103^ 


0.9 ± 0.050 


1.0 ±0.14 


0.9 ± 0.053 


0.9 ±0.13 


J 1606+4751 


54 ±8 


109.0 ± 4.4 


670 ±44 


1660 ± 166^ 


1.4 ±0.061 


0.6 ±0.13 


1.2 ± 0.050 


1.0 ±0.14 


J1618+5210 


45 ±7 


112.3 ±4.5 


294.0 ± 20 


590 ± 29.5' 


0.7 ± 0.063 


0.7 ±0.13 


0.7 ± 0.032 


0.8 ± 0.096 


J 1629+4937 


87 ± 11 


197.0 ± 7.9 


499 ± 24 


770 ± 103^ 


0.7 ± 0.052 


0.7 ±0.11 


0.6 ± 0.064 


0.5 ±0.16 


J 1637+3223 


49 ±8 


147.5 ±5.9 


438 ± 24 


660 ± 1032 


0.8 ± 0.055 


0.9 ±0.14 


0.7 ± 0.074 


0.5 ±0.19 


J1641+4209 


67 ±9 


263.7 ± 10.5 


1186 ±25 


1810 ± 181^ 


1.1 ±0.040 


1.1 ±0.12 


0.9 ± 0.050 


0.5 ±0.12 


J 1643+45 18 


40 ±7 


110.4 ±4.4 


304 ± 20 


480 ± 1032 


0.8 ± 0.061 


0.8 ±0.15 


0.7 ± 0.099 


0.5 ± 0.25 


J171 1+3047 


41 ±7 


124.8 ±5.0 


407 ± 35 


670 ± 1032 


0.9 ± 0.074 


0.9 ±0.14 


0.8 ± 0.073 


0.6 ± 0.20 



Fro dHales eTaTlfMl 
Fron lHales eTaTlfT^gTI 



TABLE 6 

Differential is'-BAND Number Counts 



K 


Galaxies 


Stars 


11.625 


1 


42 


11.875 


1 


54 


12.125 


4 


66 


12.375 


4 


77 


12.625 


8 


87 


12.875 


8 


119 


13.125 


11 


155 


13.375 


22 


167 


13.625 


28 


184 


13.875 


32 


179 


14.125 


43 


235 


14.375 


62 


259 


14.625 


85 


339 


14.875 


117 


345 


15.125 


199 


397 


15.375 


257 


477 


15.625 


381 


494 


15.875 


518 


531 


16.125 


655 


545 


16.375 


858 


693 


16.625 


1223 


713 


16.875 


1490 


751 


17.125 


2204 


597 


17.375 


2786 


496 


17.625 


3256 


475 


17.875 


3784 


443 


18.125 


4230 


376 


18.375 


4402 


274 


18.625 


4386 


237 


18.875 


3813 


145 



DRAGONS 



TABLE 6 — Continued 



K 


Galaxies 


Stars 


19.125 


2820 


107 


19.375 


1026 


31 


19.625 


210 


10 


19.875 


191 


3 


20.125 


113 


2 


20.375 


110 





20.625 


87 





20.875 


73 


1 



